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Abstract. Three sets of ladder operators in spheroconal coordinates and their respective 
actions on Lame spheroconal harmonic polynomials are presented in this article. The poly- 
nomials are common eigenfunctions of the square of the angular momentum operator and of 
the asymmetry distribution Hamiltonian for the rotations of asymmetric molecules, in the 
body-fixed frame with principal axes. The first set of operators for Lame polynomials of 
a given species and a fixed value of the square of the angular momentum raise and lower 
and lower and raise in complementary ways the quantum numbers ni and n2 counting the 
respective nodal elliptical cones. The second set of operators consisting of the cartesian 
components L^,, Ly, of the angular momentum connect pairs of the four species of poly- 
nomials of a chosen kind and angular momentum. The third set of operators, the cartesian 
components Px, Py, Pz of the linear momentum, connect pairs of the polynomials differing in 
one unit in their angular momentum and in their parities. Relationships among spheroconal 
harmonics at the levels of the three sets of operators are illustrated. 
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1 Introduction 

One of our contributions to the Symposium on Superintegrability, Exact Solvability and Special 
Functions dealt with the topic of "Symmetries and asymmetries in quantum systems confined 
by elliptical cones" , reviewing some of our results on the rotations of free asymmetric molecules 
[9, 10] and their extensions for other systems under confinement [17, 18, 19]. While the free 
systems involve Lame spheroconal harmonic polynomial eigenfunctions, the systems confined by 
elliptical cones involve quasi-periodic Lame eigenfunctions represented as infinite series, due to 
the breaking of the parity symmetry with respect the axis of the confining cone. Nevertheless, 
the methodology for the accurate and convergent evaluation of eigenvalues and eigenfunctions is 
the same. The differences reside in the boundary conditions, translating into recurrence relations 
of three terms and four terms for the expansion coefficients, respectively; such relations can be 
cast into matrix forms which are of finite size for the free systems [10] and of infinite size for 
the confined systems [17, 18, 19] leading to the corresponding polynomial and infinite series 
representations. The latter can be evaluated with sufficient accuracy using finite and large 
enough size matrices, testing for convergence along the way. 

Concerning the conditions of exact solvability requiring polynomial solutions and ladder ope- 
rators connecting them, the Lame spheroconal harmonics for the free systems satisfy the first 

*This paper is a contribution to the Special Issue "Superintegrability, Exact Solvability, and Special Functions" . 
The full collection is available at http://www.emis.de/journals/SIGMA/SESSF2012.html 
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one and in this contribution three sets of ladder operators in the purely spheroconal formulation 
are reported. It must be pointed out that one of these sets had been identified sometime ago 
[13, 14, 15, 24, 25] and our recent work [18] in its section "On Developing the Theory of Angular 
Momentum in Bases of Lame Spheroconal Harmonics" identified a second set. However, it 
must also be admitted that in these treatments the connections with the spherical and cartesian 
harmonic bases were used. 

The material in the manuscript is distributed in the following way: Section 2 reviews briefly 
the simultaneous separability of the square of the angular momentum operator and the asym- 
metry distribution Hamiltonian eigenvalue equations, the integration of the Lame differential 
equation in the respective elliptical cone coordinates, illustrating the lower eigenvalues and eigen- 
functions for their matched products making up the Lame spheroconal harmonic polynomials. 
Section 3 identifies the raising and lowering actions on the numbers ni and n2 counting the 
nodal elliptical cones in the neighboring spheroconal harmonic polynomials as their respective 
eigenvalues h^_^{ki) and h^,^{k2) are complementarity increased and decreased, which follows as 
a corollary to the tridiagonal matrix diagonalization method of solution, for chosen values of £ 
and species AB of the polynomials. Section 4 identifies the shifting actions of the Lx, Ly and Lz 
operators on the individual Lame polynomials of species A and B, respectively, and on their 
matched products of species AB for each fixed value of i. Section 5 identifies the correspond- 
ing actions of the px, Py and Pz operators on an initial harmonic polynomial iAB, raising and 
lowering the angular momentum by one unit £' = £ it 1 into final polynomials i'A'B' of oppo- 
site parity and different species. Section 6 discusses the main results of this work, including 
relationships among the spheroconal harmonics in the three levels of the identified operators, 
as well as connections with other works in the literature. 



Kramers and Ittmann [6] pioneered the quantum mechanical study of the rotations of asymmetric 
molecules using spheroconal coordinates, showing the separability of the eigenvalue equations 
for the square of the angular momentum and the complete Hamiltonian 



in the body-fixed frame with principal axes, involving the moments of inertia Ii < I2 ^ Is- The 
separated equations are of the Lame type, and its polynomial solutions are exact in principle, 
but could not be numerically implemented for higher excited states at that time. Consequently, 
the study of the rotations of asymmetric molecules developed along the perturbation theory 
route taking the symmetric prolate and oblate molecule spherical harmonic solutions as the 
non-perturbed starting point [7, 8, 16, 29, 35]. 

Patera and Winternitz used an alternative energy operator 



in their work [24]. On the other hand, Piha [27] and Valdez and Piiia [32] used the alternative 
parametrization of the original Hamiltonian of equation (1), in the form 



2 Eigenvalues and eigenfunctions of L and H* 





E = -A{L\ + rLl) 



< r < 1 



H = -QL" + -P{eiLl + e2Ll + e^Ll) 



where 
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is the average of the three inverses of the moments of inertia characterizing a spherical top, 
2 



p2 



1 /I /I 1^^ 
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measures the magnitude of the asymmetry of the molecule and 

H* = ]^{eiLl + e2Ll + e^Ll) (2) 
is the asymmetry distribution Hamiltonian, with parameters ei, ei and 63 such that 

ei > 62 > 63, 61 + 62+63 = 0, ei+62 + 6| = ^. (3) 

The inequalities follow from those for the moments of inertia, the vanishing of their sum 
reflects the zero trace of the matrix of the inverses of the moments of inertia after the term in Q 
is separated, and the coefficient in and the sums of the squares are correlated. Only one of 
the parameters 6j is independent, the set of Ii, I2, I3 being replaced by Q, P, 6j. 

Both of our works [9, 10] used the parametrizations of equations (2) and (3) with a single 
independent parameter 6j, in analogy with that in equation (1). At that time we were not aware 
of [24]. Not surprisingly, their equivalence can be established and our respective methodologies 
also overlap. In fact, our solution using the spherical canonical bases to construct the matrix 
representation of H* and determine its eigenvalues and eigenfunctions by diagonalization [10] is 
related to Sections 2 and 3 in [24]. On the other hand, our construction of the Lame eigenvalues 
and polynomials also by diagonalization of the matrices representing the recurrence relations 
for the series expansions in the Jacobi elliptic function representation of the spheroconal coor- 
dinates [10], has as well points of contact with parts of Section 6 in [24]. We have just become 
aware of some additional references [13, 15] and [14] of that period, whose titles describe the 
overlap with this contribution. Specifically, the role of the D2 group of rotations by vr around 
the cartesian coordinate finite symmetry group is explicitly recognized in [25] as the 

missing ingredient for the complete characterization of the different species of Lame spheroconal 
harmonics. 

We refer the reader to [10, 27, 32] for the details or equations in what follows. The transfor- 
mation equations between spheroconal and cartesian coordinates are 



X 



rdn (xij^i) sn (x2|A;2), y = rcn (xil^i) cn {x2\k2), 



z = rsn{xi\ kf) dn (x2 1 kj) (4) 

in terms of Jacobi elliptic functions sn [xlkf) , cn {xlkf) and dn [xlkf) with matching parameters 
kf -\- k2 = 1 [1, 20]. The elliptic functions are related as 

cn2(x|A:2) =l-sn2(x|A:2), dn'' {x\k^) = I - k^ sn^x\k^) , (5) 

and their derivatives are 

-^sn (xW) = cn dn ^ cn {xW) = - sn (x|A;^) dn (xW), 

^ dn (xl^^) = -k'^ sn {xW) cn (xl^^) • (6) 

The transformation equations (4) and equations (5) allow the identification of r as the radial 
spherical coordinate, xi a-s an elliptical cone coordinate with axis along the x-axis and X2 as an 
elliptical cone coordinate with axis along ^;-axis. Their respective domains are — tt < am(xi) < tt 
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and —11/2 < am(xj) < 7r/2 for kf > k'j. Since k\ is always associated with variable xi and ^2 
with from now on in this manuscript, we simplify the notation as {xi\kf^ — t- (Xi)> ^^^^ 
sake of space saving. 

The transformation equations (4) and the equations (6) lead in turn to the scale factors 



hr = l, h^ = h^^=h^^=r^Jl-klsn^{xi) -klsv?{X2) (7) 
and to the Laplace operator 

^2 1 9 / 2 5 \ 1 / 92 52 



V^ = — — r — +— - h 

From the latter the identification of the square of the angular momentum is immediate: 

The cartesian components of the angular momentum vector itself can be found in the Appen- 
dix of [10], and allow the construction of the explicit form of the asymmetry distribution Hamil- 
tonian 



^* = (^(ei - (ei - 62) ST?{X2)) ^ + (es - (e2 - 63) sn2(x2)) 



(9) 



Since L? commutes with H* ^ as easily tested using equation (2), both operators admit common 
eigenfunctions with the factorizable form 

^'(xi,X2) = A(xi)A(X2) (10) 

for which the eigenvalue equations of the operators of equation (8) and (9), with eigenvalues 
+ 1) and E* become separable. Indeed by using the angular scale factors of equation (7), the 
separability can be accomplished provided the geometric parameters in equations (4) and the 
dynamic parameters in equation (9) are connected by 

,2 _ 62-63 2 _ ei -e2 

ei - 63 ei - 63 

This means that the specific spheroconal coordinate system is tailored to the asymmetry distri- 
bution of the molecule. 

The next result is that each of the factors in equation (10) satisfies the Lame equation 
[21, 33, 34] 



^+l{^ + l)kfsn\xi) 



Kixi) = KKiixi) (11) 



for i = 1,2, with the respective eigenvalues 

= + ,, = !(i±iM_^. (12) 

ei - 63 61 - 63 61 - 63 61 - 63 

Their combinations 

hi + h2 = l{^ + l), eihi + e3h2 = 2E* , (13) 
provide the eigenvalues of and H* , respectively. 
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Table 1. Matching Lame polynomials according to species and parities. 



A 

A 


1 


a 


c 
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sc 


sd 


cd 


scd 


AB 


1 
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z 
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xz 


yz 


xyz 


n,. 


+ 




+ 


+ 






+ 




n. 


+ 


+ 




+ 




+ 






n. 


+ 


+ 


+ 




+ 









Next, the solutions of the Lame equation (11) are also reviewed. The derivative equations (6) 
for the elliptic functions, represented by their initial letters allow the identification of 1, s, c, 
d, dc, ds, sc, dcs as removable singularity factors in the differential equation. The dependence 
of both operators and H* on the squares of L^, Ly and Lz guarantee that their common 
eigenfunctions, equation (10), have definite inversion parities under x — ?• —x, y — t- —y, z — )• —z. 
Correspondingly, the Lame functions get classified into two kinds, each with four species sharing 
the same number of even or odd factors in the removable singularity factors. The Lame functions 
in the respective elliptic cone coordinates xi and X2 in equation (10) must be matched according 
to their singularity factors A(xi) and B{x2), and the respective parities from equations (4), as 
illustrated in Table 1. 

There are two kinds of Lame spheroconal harmonic functions distinguished by their over- 
all parity (— )^ = IlajIIyn^. For each value of i each kind may contain four different species 
characterized by the respective singularity factors. 

The solutions of the Lame differential equation (11) have the general form 

' max 

AHx^) = Aix^)^afsn'■^iXi). (14) 

s=0 

as a series of even-power of the elliptical sn{xi) function, where the expansion coefficients 
satisfy three-term recurrence relations, which the interested reader may see in equations (32)- 
(39) of [10], for the eight different species A. For the lower i eigenstates the eigenvalues h^[kf^ 
and ratios of the expansion coefficients a^/a^ can be evaluated in a straightforward way from 
equation (11), and the linear independence of the powers of sn{xi)- The corresponding results 
are illustrated in Table 2 for the individual Lame polynomials and in Table 3 for the matched 
spheroconal harmonics. 

For higher values of £ it is more practical to cast the recurrence relations into eigenvalue 
problem matrix forms and to obtain from their diagonalizations the respective eigenvalues and 
expansion coefficients. The sizes of the matrices are as follows for i even: 

= - + 1 and N"" = N""^ = N'"^ = - 

max 2 t^ii--! i max max max 2 ' 

and for ^ odd: 

/ _i_ 1 f — ^ 

AT* = TV = N'^ = and N'''"^ - 

' max ' max ' max o cii-iu. j. r ^^^-^ 



In both cases the total dimensions are 21 + 1, corresponding to the number of independent 
spheroconal harmonic polynomials for each value of i. For each species there are -/V^ax different 
eigenvalues which can be ordered according to their increasing values 



< < < ■ ■ ■ < hf < ■ ■ ■ < h^A .<h^A , (15) 

* max * max 
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or their decreasing values 

> >--->hf>--->hi>h^>hf. (16) 

Obviously, the expansion coefficients in equation (14) depend on the specific eigenvalues as{hf), 
as{h^)- The first of equations (13) suggests the matching of the smaller eigenvalues of one 
species with the larger ones of the companion species for the construction of the complete set of 
spheroconal harmonics iABij. 

Let us also analyze the number of nodes in the respective factors of Lame functions in equa- 
tion (14), and n^. For the singularity removing factors their specific values are as follows: 

= n'' = 0, n^^ = or 1, = 1, 

n'^" = or 1, n'^' = 1, rf' = 1 or 2, r/^' = 1 or 2. (17) 

The polynomial of degree 2N^g^^ in sn(x) has at the most N^.^-^ zeros sn^(xp), appearing by 
pairs with opposite signs it sn(xp) with the possible numbers 

n?' = 0,2,4,...,2(iV^ax-l) (18) 

counting the symmetric sheets of the nodal elliptical cones. 

Therefore, the total number of nodes for A.^^{xi) Lame polynomial is 

m =n'^ + nP = 0, 1,2,...,^ (19) 

and likewise for its matching companion A^^(x2) 

n2 = n^ + nP = 0, 1,2,...,£ (20) 

Correspondingly, the spheroconal harmonic functions from the product of the pairs of matched 
Lame polynomials exhibit the possible numbers 

n^ + n^ = 0,1,2,3 (21) 

counting the shared nodal cartesian planes, as limiting cases of nodal elliptical cones, for the 
respective species 

AB = 1\ X, y, z; xy, xz, yz; xyz. 

In turn, 

ni + n2=e (22) 

counts the total number of nodes characterizing any harmonic function of order i. 

The ordering indices i, j for the eigenvalues in equations (15) and (16) can be replaced with 
the quantum numbers ni and 77-2. This is completely justified by recognizing that the eigenvalues 
with matching quantum numbers, according to equation (22), are also matched to satisfy 

/^If + C = + 1) (23) 

as required by equation (13). 

Table 2 illustrates the eigenvalues hf^ of the Lame functions for the successive values of 
^ = 0, 1, 2, . . . , and the respective species with the singularity removing factors. The eigenvalues 
and their quantum numbers ni = 0,1,2,... are listed in increasing order. The associated 
polynomial factors in the eigenfunctions, equation (14), are also included. 

Table 3 illustrates the matching of Lame eigenvalues and polynomials to obtain the sphero- 
conal harmonics iABnin2 and their eigenenergies 2E* from equation (13), written in terms of 
the asymmetry distribution parameter Cj. The readers may check the ordering of the individual 
Lame eigenvalues, correlated with the increasing of the number of nodes in the polynomials in 
Table 2. The matching according to equations (22) and (23), and the ordering of increasing 
energies can be appreciated in Table 3. 
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Table 2. Lame eigenvalues and polynomials. 
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Polynomial 


£ = 
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riQ — 
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1 
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r rl 


7 9 


1 






c 


hi = 


1 


1 






s 


hi = 


1 + ^1 


1 


i = 


2 


dc 


hf = 




1 






ds 


hf = 


-I , A 1 

= 1 + 4fcf 


1 






cs 


hf = 


--i + kf 


1 






1 


hi = 


2(l + fc2)-2^1-fc2(l-fc2) 
V 


1- ^sn2(x.) 






1 


hi = 


2{l + q) + 2^l-kl{l-k^) 


1 - ^sv?{xi) 


£ = 


3 


dcs 


hf' 


= 4(1 + fc?) 


1 






d 


K = 


(5fc2 + 2)-2^4A:4-A;2 + l 


1 + ^^sn2(xi) 






d 


hi = 


(5fc2 + 2) + 2^4fc,4-A;2 + l 


l + ^sn2(x.) 






c 


hi = 


(5 + 2fc2)_2^4_^2(i_fc2) 


l + i^sn2(xO 






c 


h% = 


(5 + 2fef) + 2^4-A;f(l-A;2) 


1+ i^sn2(xi) 






s 


hi = 


5(l + fcf)-2^4fc4-7fc2 + 4 


1+ 'sn^(xi) 






s 


hi = 


5(l + fc2) + 2y/4fc4-7fc2 + 4 


1 + I±^sn2(x,) 


f) 

I = 


4 


1 


h — 


20(1 + k^)h^ 

+ [D4(,i + ; + -^Uoftj JAi 

— O^UfVj ^1 "T A/j j — U 


l+^sn2(x,) + ^sn4(x,) 

± ne rooxb oi xne polynomial 

are /iq, /i^, /14 






dc 


/if = 


--b{l + kj)-2^A + k-i + Akl 


1 - — ^sn^(xi) 






dc 


hf = 


= 5(l + fc2) + 2^4 + A;2 + 4A;| 


1 '^3-a+^?) 3n2(xO 






ds 


/if = 


= 5(l + 2/cf)-2y/4-9fc2(l-/c|) 


1 '^"-r^') -2(x.) 






ds 


hf-- 


= 5(l + 2A;2) + 2^4-9fc2(l-A;2) 


1 '^3--(l+4.f)^^2(^^) 






cs 


hf = 


= 5(2 + kf) - 2^9 - 9/^2 + 4/c4 








cs 


hf = 


= 5(2 + kf) + 2^9 - 9A;| + 4/c4 


1 - sn2(x.) 



3 Raising and lowering of quantum numbers rii and n2 
for spheroconal harmonic polynomials of a chosen kind 
and species iABnin2 

The spheroconal harmonic polynomials with their labels for angular momentum i, species AB 
and number of nodes rii and 112 can be written as the products of the pairs of matching Lame 
polynomials 

<n^.(xi,X2;/.lf,0 = Alf (xi,0a1?(x2,0> (24) 

^ max 

Alf (Xi, O = ^(Xi)ao^(O E afiKi) sn'^{xi) (25) 

s=0 
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Table 3. Eigenenergies 2E* in units of fi^ for Lame spheroconal harmonics £ABnin2- 
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scd 
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1 
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2 
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1 


—362 — 3y 5 — 462 
_3eo - 3,/5 — 4e? 





d 


g 


2 


1 


-361 + 3^/5 - 46? 




c 


c 


3 





-362 + 3^5-46^ 






d 
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_3p„ _|_ 3. /5 _ 4p2 


A 

1 = 4: 


■i 

1 


■i 

1 


U 


4 


eih^[k{) + 63/i|(fc|) 




dc 


sc 


1 


3 


063 — y 21 + 0U63 




ds 


sd 


1 


3 


562 - + 606^ 




cs 
1 


ca 
1 


o 
z 

2 


o 
z 

2 


O61 — Y Zi + DU6]^ 

6i/i^ (A;?) + 63/12(^1) 




dc 


sc 


3 


1 


563 + V2I + 606^ 




ds 


sd 


3 


1 


562 + + 606^ 




cs 
1 


cd 
1 


4 
4 






561 + \/21 + 606^ 
eihi{kl) + e^hlikl) 



factorizing the zero power coefficients (h^ril) using the ratio of coefficients for the other 
powers 



A 



= ^TO^Y (26) 



and similarly for its companion with A ^ B and ni — )• 77.2- As established in the previous section, 
equations (22) and (23) define the matching of quantum numbers and eigenvalues, respectively. 

For chosen values of £, A and B, the removed singularity nodes are fixed in numbers 
and , according to equations (17)-(21). Then according to equations (18)-(20) the number 
of nodes n^, ni and n2 change in steps of 2. the inequalities in equations (15) and (16) for the 
increasing and decreasing eigenvalues become 

"n^ < "'n'4+2 < • • • < "'ni-2 < < "ni+2 < " " " < 'i£_„s_2 < '^i-nB^ l^' J 



and 



> '^£-n'4-2 > " " " > ^2+2 > > K2-2 > ■ > Kb +2 > Kb- (28) 



If we start from the state lABnin2 of equation (24) depending on the eigenvalues appearing 
in the middle of equations (27) and (28), the state involving their neighbors to the right have 
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raised ni by two units and lowered n2 by two units with the corresponding change in the 
eigenvalues which remain satisfying matching conditions of equations (22) and (23); the same 
changes operate in equations (25) and (26) for the eigenfunctions via the eigenvalue dependence 
of the expansion coefficients. Similarly, if we had moved to the neighbors on the left: ni — )■ ni — 2 
and 712 — 7- ni + 2 with the complementary lowering and raising by two, with their consequences 
in the changes in the matching eigenvalues and the other state of the same species. Obviously, 
the lowering and raising actions can be extended to the right and left ends of equations (27) 
and (28) covering all the N^^^ = N^^^ states of the lAB species. Examples of these connections 
can be found in Tables 2 and 3. 

As a prelude to the following two sections we include a couple of purely spheroconal coordinate 
tools. The first one illustrates the successive species changing actions of the derivative operator 
on the singularity removing factors and even-power polynomials appearing in the process: 





a 






dxi 


9x1 


1 








s 


cd 


{-kfc^ - d^)s 


c 


—sd 


{-kfs^ + d2)c 


d 


-kfsc 


-kf{c^ - S2)d 


cd 


{-kfc^ - d^)s 
{-kfs^ + c^)c 


(-fc2c2 - d2 + 2kfs^)cd 


sd 


{kfs^ -d^- 2kfc^)sd 


sc 


(s2 - ^)d 


{-kfs^ + kfc"^ + 2kfd'^)sc 


scd 


-k'^c^s'^ - d'^s'^ + d^c^ 


2[kf{s^ - C2) + {kfs^ - d2) _ (^2^2 + d^^^scd 


s\ c\ d^ 


2scd{l,-l,-kf) 





The first derivative causes changes of species with one factor to species with the other two 
factors and viceversa. The species 1 is represented by the top entry, and also by the bottom 
one; and species dcs and 1 are also connected by their first derivatives as illustrated by the two 
bottom entries. The second derivative returns each species to the original one. The even-power 
polynomials can be reduced to polynomials in sn^lx) via equation (5). 

The second tool is connected with equation (25), expressing the set of Lame polynomials £Ani 
as a linear combination of the A(x) sn2^(x) basis. Both bases have the same dimension N^^^^ + l, 
and the members of the latter can be expressed as linear combinations of Lame polynomials. 
The transformation coefficients af^^ are the elements of the inverse matrix of a^{h^). 



4 Angular momentum cartesian components connecting pairs 
of polynomials of the four different species with a common 
eigenvalue i 

The explicit expressions for the cartesian components of the angular momentum operators in 
spheroconal coordinates are borrowed as equations (A15)-(A17) in [10]: 



h 



dn(xi)cn(x2)dn(x2 



_d_ 
d 



+ kl sn(xi)cn(xi)sn(x2) 



d 
d 



-cn(xi) sn(x2) dn(x2)7^ h sn(xi) dn(xi) cn(x2) 

dxi OX2 



-kl sn(xi) sn(x2) cn(x2) 



d 
dxi 



cn(xi)dn(xi)dn(x2 



d 
9X2 



(29) 
(30) 
(31) 
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Table 4. Eigcnfunctions with £ = 1 and its resulting application of each cartesian angular momentum 
operator. 





Lx 


Ly 


Lz 


A^Af 
AfAg 
AfA^ 




AfA^ 
-AfAg 




AoAf 


AfAe 

-A^Af 




This section evaluates the actions of these operators on the lower iABnin2 spheroconal 
harmonic polynomials of the same kind for a fixed value of £, connecting with different species 
eA'B'n[n'2. 

The calculations are straightforward requiring care, labor and some guidelines. The deriva- 
tives d/dxi act only on the first or second factor in equation (24) for i = 1 or 2. The individual 
changes of species were already described in the first tool of the prelude, and for an initial mat- 
ching species AB the final ones A'B' follow. Notice the appearance of the square of the angular 
scale factor in the three operators in equations (29)-(31). This requires that the application of 
the operators inside the brackets to the spheroconal harmonic polynomials provides the A'B' 
factor and also the /i^ /r^ factor to compensate the one in the previous line. This requirement is 
indeed satisfied. The remaining factor can be reduced to a polynomial in even-powers of sn^(xi) 
and sn^(x2), up to A^^ax = -^m^x) respectively. The remaining task is to identify the harmonic 
linear superposition of the polynomials, making use of the linear independence and complete- 
ness of both bases in the second tool of the prelude. Some of our results on the actions of Lx, 
Ly, Lz on the non-normalized spheroconal harmonic polynomials iABnin2 for ^ = 0, 1, 2, 3 are 
illustrated and described next. The application of any of the operators Li on the AqAq = 1 leads 
to the same single eigenstate with its eigenvalue zero. 

Table 4 exhibits the result of the corresponding applications on the three spheroconal harmo- 
nics with £ = 1. The contents in the table should not be surprising upon recognizing the species x, 
y, z of the successive monomials, and their equivalence with the respective cartesian harmonics. 

Table 5 for the i = 2 polynomials contains familiar information for the unique monomial 
species xy, xz and yz transforming among themselves under the respective rotations, with 
coefficients one in the corresponding entries. The two AqA2 and AgAg binomial companions get 
transformed into the monomials of the other species as illustrated by the top and bottom rows 
in the table, with coefficients involving the eigenvalues hl^_^{kf) and /i.^2(^2) matching 
binomials; the respective species of the monomials are complementary to that of the operators, 
as the reader can ascertain by moving along the top and bottom rows in the successive columns. 
In turn, the respective operators acting on the monomials lead to the linear superpositions of 
the two A^^A,^2 binomials with more elaborate coefficients including the common one 

^1 -2-ai{hl{kl)) + -4{ht{kl)) 2 

3 = Cli (by changing d^s), 3 = - (l + C], 3) , (33) 

^1 ^_^2 al{hl{kl)) + -al{hl{kl)) 

af(/i|(fc2))^c(;,c(^2))_^c(;,c(fc2))^c(;,c(fc2))- I ) 

The spheroconal harmonics with I = 2> come in a single monomial of species [(ics][,scd] = 
[xyz] and three pairs of binomials of species x, y and z. The monomial appears in the middle 
row of Table 6, and its companions of the other species in the rows below and above. The 
monomial transformed into the x, y, z species in the successive columns by the respective 
angular momentum components; the transformed states are linear combinations of the pair of 
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Table 5. Eigenfunctions with £ — 2 and its resulting application of each cartesian angular momentum 
operator. The common factor in the table has the value — [/ig(fcj) — hiikf)^ 

Lx Ly Lz 



Af Af 
Af Af 



Af Af 

-Af Af^ 

{2-hi{kl))KiKi 
{2-hi{kl))KiAi 
-hi{kl)KfKf 



{hi{kl)-hi{kl))KfKf 



2C^ 



-KM 



Af Af^ 
{hi{kl)-hi{kl))KfKt 



/i^(A;f)Af Af 
{2-hi{kl))KiKl 
{2-hi{kl))KlKi 

-Af Af 
hl{kl)KfKf 



Table 6. Eigenfunctions with ^ = 3 and its resulting application of each cartesian angular momentum 
operator. The element ^ represents the k coefficient of the linear combination at the i row and j 
column of the table. Some examples are represented in the equations (32)-(34). 



Lx 



L,. 



A^A| 
AfA^ 
AfA^ 

\dcsj^scd 

A^Af 
A§Ag 
A|A^ 



2af {h^ikf)) Af ^Af 
AfAf + C|,A|A; 



cd 



^1 ACAC 





+ C^iiAiAg 



^^i,iA^A| + Cl.AlAl 
2ai {h^{kl)) Af '^Af 
CiiAfA^ + C^^AfA^ 
C^^iA^A^ + C^^AiAg 



AfA^ + C2 A|A^ 



ikl))Ai 



cs ^scd 



A^Af 



^1 ACAC 



+ C'4_2A§A[j 



Cl,AlAi + Cl,AiAd, 



/io(^2))Af'Af^ 



Cl 2A0AI 



A^Af 



C'l.sAfAg + (7^3 



CIsAqAI 



+ CLAiAl 



-2ai (h^ikj)) Af'Af^ 

C} 3A^A2 + CjgAqAn 



4,3AfAi 
CisAoAl 



+ ch 



A'^ A'^ 



+ C2 

2ai [h^^iki)) Ar Af 



Ad AS 

dcs \scd 



companion binomials with different number of nodes ni and n2- The inverse transformations 
of the successive binomials leading back to the single monomial are identified in the row and 
column positions (1, 1), (2, 2), (3, 3), (5, 1), (6, 2) and (7, 3). The individual binomial spheroconal 
harmonics are transformed also into linear combinations of the pairs of binomials of species 
complementary to the component of the angular momentum, as an inspection of the remaining 
twelve positions in the table show. The coefficients involved become more numerous, taking 
into account the increasing number of the states and their combinations. 

We complement this section with the following diagrams in Fig. 1 sketching the connections 
among the four species, for the two kinds of even and odd i of spheroconal harmonics, provided 
by the angular momentum components. The interested readers may compare them with their 
counterparts in [24]. The coefficients in the entries of Tables 4-6 are the counterpart of the 
raising and lowering operator coefficients for the spherical harmonics. 




< ► 3^ 



xyz 



Figure 1. Action of the angular momentum operators on each of the spheroconal harmonics. 
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5 Linear momentum cartesian components raising and lowering 
the angular momentum by one unit between polynomials 
of opposite parities 



The cartesian components of the hnear momentum operator in spheroconal coordinates are 

ihr 



Px = dn(xi)sn(x2)Pr 
Py = cn{xi) cn{x2)Pr 
Pz = sn(xi)dn(x2)Pr 



hi 
ihr 



hi, 
ihr 



h sn(xi)cn(xi)sn(x2 
-sn(xi)dn(xi)cn(x2) 
cn(xi)dn(xi)dn(x2 



+ dn(xi)cn(x2)dn(x2^ 



d 

d .2 



dxi 



cn(xi)sn(x2)dn(x2) 
H sn(xi)sn(x2)cn(x2) 



9X2. 
d 
9X2 
d 
9x2 



,(35) 
(36) 
.(37) 



Notice that the angular factors of pr are simply the projections of the unit radial vector along 
the respective axes, equation (4). The angular factor of the derivatives inside the brackets in 
equations (29)-(31) and (35)-(37) for the successive components are the same but in exchanged 
positions, and different signs in some cases. The reason for this is found in the decomposition 
of the p vector into its radial and transverse components 



p = rpr — r X 



(38) 



already used in [18]. Correspondingly, also notice the presence of the factor r/h^ in equa- 
tions (35)-(37). The actions of these on a spheroconal harmonic functions £ABnin2 lead to 
functions of the other kind i'A'B'n\n2 with i' = £±1, species A'B' with the other parity. Con- 
sequently, the method and some of the steps and results of the previous section are useful or have 
their counterparts in the application of the operators of equations (35)-(37) on the spheroconal 
harmonic polynomials iABnin2. 

At the level of each individual spheroconal coordinate, the multiplications by d, c, and s 
of the first terms in the operators in equations (35)-(37) with the eight species of singularity 
removing factors ^(x) lead to the following result: 





1 


s 


c 


d 


cd 


sd 


sc 


scd 


d 


d 


ds 


dc 




cd' 


sd' 


scd 


scd' 


c 


c 


cs 




cd 


c^d 


scd 


sc' 


sc'd 


s 


s 




sc 


sd 


scd 


s'd 


s'c 


s'cd 



with the common result that every original A{x) is changed in kind and species. The effects 
are qualitatively similar and complementary to those of the first derivative in the table of the 
prelude at the end of Section 3. In fact, the union of both tables provides the sets of four species 
associated with the other kind of Lame functions. 

What has been said for A{xi) also holds for the companion B{x2), and the resulting ^'(xi) 
and B'(x2) can be properly matched. At the level of multiplication in the first terms of 
equations (35)-(37), their respective x, y and z character has already been recognized. As 
a consequence of the exchange of angular factors in the derivative terms in equations (29)-(31) 
and (35)-(37), the kinds and species maintained in the previous section are changed in the 
present one. The net result is that the angular derivative terms in equations (35)-(37) acting 
on the individual Lame functions iAni and £Bn2 transform them into the same kinds and 
species A'(xi), -B'(x2) obtained by multiplication in the respective first terms. 
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We proceed to analyze the effects of the operators of equations (35)-(37) on the successive 
Lame spheroconal functions. For i = 0, Aq(xi)Aq(x2) = 1 is annihilated by the derivative terms, 
and transformed into 

dn(xi) sn(x2) = Ao(xi)AJ(x2) = ^, cn(xi) cn(x2) = A5(xi)Ag(x2) = ^, 

sn(xi)dn(x2) = A^i(xi)A^(x2) = ^, 

respectively, by the multiplication operations. The identifications with the cartesian harmonics 
is also made. Their angular momentum of £ = 1 and negative parity are also recognized in 
contrast with the i = and positive parity of the original state. 

Next, we illustrate the net result of the applying the derivative terms on the i = 1 states, 
reducing them by factoring the square of the scale factor, obtaining the following table of mono- 
mials and binomials: 

Table 7. Actions of angular derivative terms in px, Py, Pz on i — 1 spheroconal harmonics. 





Px 


Py 


Pz 


A^'^A^'^ 
Aj^A^'^' 


l-dn2(xi)sn2(x2) 


A_2dc^2sc 

1 -cn2(xi)cn2(x2) 

\ 2cs A 2cd 
^^2 ^^0 


^2ds J\^2sd 

\ 2cs A 2cd 
^^2 ^^0 

l-sn2(xi)dn2(x2) 



The off-diagonal entries are identified as A|'^'^(xi)A^'''^(x2) = xy, A|'^'*(xi)Af*'^(x2) = xz and 
A|'^^(xi)Aq'^'^(X2) = obtained via two different operators acting on two different initial states. 
They could also be obtained via the multiplication with the factors in the first terms of equa- 
tions (35)-(36). The identification of their angular momentum £ = 2, and positive parity and 
nodal quantum numbers can be made from Table 3. 

The diagonal entries deserve additional analysis before the lowering effect of the operators 
is identified, as well as the raising effects leading to the remaining two i = 2 spheroconal 
harmonics. The three entries have been obtained via three different operators acting on three 
different initial functions. Of course, the "ones" correspond to the monopolar harmonic to which 
we go back by the lowering effect on the dipole harmonics. On the other hand, the other terms 
are identified as the squares of the projections of the radial unit vector along the coordinate 
axes. Their sum is "one" coinciding again with the lowered £ = state. Of the three terms, 
only two are linearly independent, since their sum is fixed. This suggests that the three terms 
can be combined in two linearly independent combinations to be identified with the missing 
Aq^(Xi)A2"'(x2) and A2^(xi)Aq^(x2)! £ = 2 spheroconal harmonics. 

The analysis can be made more transparent by recognizing that x^, y'^ and z'^ are equivalent 
to them upon multiplication by r^. Let us consider the linear combination 

^{x,y,z) = ax"^ + I3y^ + jz'^. (39) 

It is advantageous to choose it as harmonic 

V^^* = 2a + 2/3 + 27 = 0, (40) 

because in such a case we are sure that their connection with the spheroconal harmonics is 
linear. In fact, from equation (40) we get 

/3 = -a -7, 

and rewrite equation (39) in the separable form in spheroconal coordinates 



^{r,Xi,X2) = [adn^(xi)cn^(X2) - (a + 7) cn^(xi) cn^(X2) + 7 sn^(Xi) dn^(X2)] • (41) 
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The question to be answered is which values of the coefficients a and 7 combining the diagonal 
angular entries lead to the £ = 2 spheroconal harmonics of species 1 1. 

The answers are obtained by rewriting equation (41) as polynomials in even-powers of sn(xi) 
and sn(x2) and comparing them with the target polynomials from Table 2. Their explicit 
forms are 

ao2 = -\+ \hi{ki) - \hi{ki), 102 = -\- \hm) + \hm), 

and 

«2o = "2 + \hm) - \hi{ki), = "2 - \hm) + \hm, 

defining the proportions of the presence of the respective spheroconal harmonics in the diagonal 
entries in the Table 7. We remind the reader that the analysis has been made for non-normalized 
polynomials. In any case, the coefficients above are the counterpart of the familiar spherical 
Clebsch-Gordan coefficients. Work on the raising and lowering from £ = 2,3,... is in process. 

6 Discussion 

Some of our works on the rotations of asymmetric molecules and the Lame spheroconal harmonics 
have taken advantage of their connections with their spherical and cartesian counterparts [9, 18], 
and others have been developed strictly within the spheroconal formalism [10, 17, 19, 18]. The 
analysis and identification of the three sets of ladder operators in this article has been purpose- 
fully implemented in the same formalism, avoiding the spherical and cartesian crutches. 

The connections among Lame spheroconal polynomials of a given kind and species, and 
different numbers of nodes have been explicitly recognized in Section 3. They share the same 
form and degree, differing in their number of nodes n^; in contrast with the classical polynomials 
having a number of nodes equal to their degree. For the spheroconal harmonics constructed by 
the product of Lame polynomials with matching number of nodes and eigenvalues, equations (22) 
and (23), their complete set is connected by their complementary raising of ni and lowering of 712. 
The inversion of equation (25) provides the basis of the powers of sv?{x) as a linear combination 
of the set of corresponding Lame polynomials, as the counterpart of the known relationships in 
the classical polynomials. 

Our other contribution in this volume of SIGMA [11] contains another example of polynomial 
eigenfunctions, for the Hydrogen atom confined by a dihedral angle in prolate spheroidal coor- 
dinates, evaluated by diagonalization of tridiagonal finite size matrices, and their connections 
by complementary changes in the numbers of prolate spheroidal and hyperboloidal nodes. 

The shifting actions of the Lx, Ly, operators among the spheroconal harmonics of the four 
species for a chosen value i, as illustrated in Tables 5 and 6, allow the generation of the corre- 
sponding coefficients in the purely spheroconal formalism. They are the counterpart of the famil- 
iar raising and lowering coefficients C±{i, m) for the spherical harmonics, and even closer for the 
spherical harmonics with definite parities. The relationships so far are for the non-normalized 
coefficients, but work is under way to identify and incorporate the normalization coefficients. 

The i raising and lowering actions of px, Py and pz among spheroconal harmonics of neighbo- 
ring kinds and species has been illustrated for the lower states, but the systematics had already 
been recognized in [18]. While writing the latter, several works on related topics [12, 30, 31] 
allowed us to recognize the complete radial and angular raising and lowering actions for the free 
particle in spherical coordinates. Now we can also extend it to spheroconal coordinates. 

Piha identified four step A£ = 4 ladder operators for spheroconal harmonics of the same 
species for the most asymmetric molecules [26]. We are investigating their possible connection 
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with the application of the pi operators four times in succession leading to the same states. In 
such a case the results are valid for any asymmetry distribution. 

The familiarity with cartesian, spherical and spheroconal harmonics is decreasing when writ- 
ten in this order. The works [9, 24] used the spherical basis as the tool to evaluate the spheroconal 
harmonics. In the process the interbasis expansion coefficients are generated for the original 
problem, and automatically also for the inverse problem. The specific example at the end of 
Section 5 shows the spheroconal quadrupole composition of the cartesian one in equation (39). 

As a complement to the comment in the third paragraph of this section, we include the 
following related references with their respective comments. The contents in [2] for the free 
Hydrogen atom in configuration space overlap with ours in [11], with emphasis on the rarely 
treated spheroidal sets, and in the present contribution, on the properties of spheroelliptic or- 
bitals, which have been so far practically ignored. The related work [3] also deals with the free 
Hydrogen atom in momentum space focusing on the hyperspherical harmonic elliptic Sturmian 
basis set in S"^. The earlier work [5] deals with the free quantum motion on the three-dimensional 
sphere ellipso-cyllindrical coordinates, distinguishing between prolate and oblate elliptic co- 
ordinates. On the other hand the more recent work [4] deals with the use of the hyperspherical 
elliptic coordinates and harmonic basis sets in the fifth dimensional sphere for the study 
of the three-body problem in quantum mechanics. Its classical astronomical counterpart [28] 
uses the relative radial coordinate of the first two particles and the radial coordinate of the 
third particle relative to the center of mass of the first two, instead of the hyper-radius and 
one of the five angular coordinates in S^. It is also pertinent to mention Sasaki's contribution 
in the Symposium, "Exactly solvable quantum mechanics and infinite families of multi-indexed 
orthogonal polynomials", and to cite [22, 23], involving Hamiltonians with tridiagonal matrix 
representations. The raising and lowering operators in our Section 3 have their counterparts for 
the respective eigenfunctions in the references of this paragraph. 
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